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Thermal  radiation  properties  of  turbulent  premixed  flames  have  received  little  attention  in  the  past 
perhaps  because  of  the  lower  radiative  heat  loss  compared  with  that  for  non-premixed  flames.  However, 
the  high-temperature  sensitivity  of  NO  kinetics  and  the  importance  of  radiation  in  near-limit  laminar 
premixed  flames  provide  fundamental  reasons  for  studies  of  radiation  properties  of  turbulent  premixed 
flames.  Reduced  cooling  airflows  in  lean  premixed  combustors,  miniaturization  of  combustors,  and  the 
possible  use  of  radiation  sensors  in  combustion  control  schemes  are  some  of  the  practical  reasons  for 
studying  radiation  heat  transfer  in  these  flames.  Motivated  by  this,  we  report  the  first  (to  our  knowledge) 
study  of  spectral  radiation  properties  of  turbulent  premixed  flames.  Measurements  of  mean,  root  mean 
square  (rms)  and  probability  density  functions  (PDFs)  of  spectral  radiation  intensities  leaving  diametric 
paths  at  five  heights  in  two  turbulent  lean  premised  methane/air  jet  flames  stabilized  using  small  H2/air 
pilot  flames  in  a  coflow  of  air  were  completed.  Measurements  of  spectral  radiation  intensities  leaving  three 
laminar  flames  were  also  completed.  These  data  were  used  to  evaluate  narrowband  radiation  calculations 
independent  of  the  treatment  of  turbulent  fluctuations.  Stochastic  spatial  series  analysis  was  used  to  esti¬ 
mate  instantaneous  distributions  of  temperature.  The  analysis  requires  the  specification  of  mean  and  rms 
temperature  distributions,  integral  length  scale  distributions,  and  an  assumption  of  exponential  spatial 
correlation  function.  We  specified  the  mean  and  rms  temperature  distributions  measured  by  calibrated 
narrowband  thin  filament  pyrometry.  A  simple  flame  and  mixing  model  was  used  to  relate  the  concentra¬ 
tions  of  C02  and  H20  to  the  temperature.  We  used  scalar  spatial  series  in  conjunction  with  a  radiation 
model  to  calculate  the  mean,  rms,  and  PDFs  of  spectral  radiation  intensities.  Overall,  the  model  predictions 
are  in  reasonable  agreement  with  the  data.  The  only  improvement  needed  is  in  the  area  of  capturing 
correlated  occurrences  of  high  temperatures  along  the  radiation  path. 


Introduction 

Radiation  properties  of  non-premixed  turbulent 
flames  have  been  widely  studied  both  experimentally 
and  theoretically  [1—6].  Non-premixed  flames  radi¬ 
ate  between  10%  and  60%  of  their  chemical  energy 
release  depending  on  fuel  type.  The  radiation  pre¬ 
dominantly  involves  emission  from  soot,  CO2,  and 
HaO  [1].  The  effects  of  turbulence  on  radiation  in 
non-premixed  flames  have  been  well  documented 
[1-6]. 

Radiation  properties  of  turbulent  premixed  flames 
have  not  received  much  attention  except  for  early 
computations  in  Ref.  [7],  where  the  effects  of  tur¬ 
bulence  on  radiation  properties  of  a  two-dimen¬ 
sional  premixed  methane/ air  flame  in  a  furnace  were 
analyzed.  However,  experimental  data  were  not 
available  for  model  validation.  The  lack  of  attention 
to  turbulent  premixed  flame  radiation  in  the  past 
could  be  because  of  the  relatively  small  radiative 
heat  loss,  at  least  for  atmospheric  pressure  condi¬ 
tions. 


Effects  of  radiation  heat  loss  on  the  lean  flamma¬ 
bility  limit  and  on  the  propagation  speeds  of  near¬ 
limit  laminar  premixed  flames  have  been  long  rec¬ 
ognized  (see  Refs.  [8,9]  and  their  citations).  These 
papers  contain  interesting  findings  regarding  radia¬ 
tion-induced  bifurcations  [8]  and  effects  of  self-ab¬ 
sorption  of  radiation  [9],  Premixed  flames  are  more 
prone  to  instability  and  their  radiation  signals  may 
be  used  in  a  control  strategy.  Motivated  by  this,  the 
present  work  is  devoted  to  measurements  and  time 
and  space  series  analysis  of  spectral  radiation  prop¬ 
erties  of  turbulent  premixed  flames. 

The  effects  of  turbulent  fluctuations  on  radiation 
result  from  the  nonlinear  dependence  of  the  Planck 
function  and  the  absorption  coefficient  on  species 
concentrations  and  temperature.  The  control  strat¬ 
egy  discussed  above  must  account  for  the  radiation 
fluctuations. 

The  present  work  consisted  of  the  following: 

1.  Mean,  root  mean  square  (rms),  and  probability 
density  function  (PDF)  of  the  spectral  radiation 
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TURBULENT  COMBUSTION — Premised  Combustion  Measurements 


intensities  leaving  diametric  paths  in  a  turbulent 
lean  premixed  methane/air  jet  flame  were  mea¬ 
sured. 

2.  The  temperature  statistics  in  the  flame  were  mea¬ 
sured  using  narrowband  thin  filament  pyrometry 
(TFP)  calibrated  using  measurements  in  laminar 
flat  flames. 

3.  A  stochastic  time  and  space  series  analysis  was 
constructed  following  Ref.  [3]  to  predict  the 
flame  radiation  statistics.  The  temperature  mea¬ 
surements  and  a  combustion/  mixing  model  were 
utilized  to  estimate  the  species  concentration 
time  series  in  conjunction  with  local,  mean  equiv¬ 
alence  ratio  measurements.  The  radiation  model 
was  independently  verified  using  measurements 
in  two  additional  laminar  flames. 


Background 

Axisymmetric  Premixed  Turbulent  Flames 

Axisymmetric  pilot-stabilized  Bunsen-type  turbu¬ 
lent  premixed  flames  identical  to  those  previously 
used  in  this  laboratory  [10]  were  considered.  Many 
others  [11-18]  have  used  this  flame  configuration. 
Many  of  these  studies  involved  temperature  and  ve¬ 
locity  measurements  [11-14,16],  Some  studies  have 
involved  visualization  of  flame  surfaces  and  integral 
length  scales  of  wrinkling  [15].  To  our  knowledge, 
work  has  not  been  reported  on  the  radiation  char¬ 
acteristics  of  axisymmetric  turbulent  premixed  jet 
flames. 

Thin  Filament  Pyrometry 

Accurate  measurements  of  temperature  time  se¬ 
ries  are  essential  for  evaluation  of  turbulent  radiation 
properties.  Pitts  [19]  has  shown  that  TFP,  originally 
introduced  by  Vilimpoc  and  Goss  [20],  has  the  po¬ 
tential  to  provide  5  to  10  K  precision  and  accuracy 
for  flame  temperature  measurements.  This  tech¬ 
nique  can  provide  a  1  ms  time  response  with  120  fim 
resolution  and  a  high  signal-to-noise  ratio  above 
1000  K. 

Temperature  profiles  in  diffusion  flames,  partially 
premixed  flames,  and  turbulent  premixed  flames 
have  been  measured  using  broadband  TFP  [21-23], 
With  a  liquid  nitrogen  cooled  InSb  detector  oper¬ 
ating  over  a  range  of  wavelengths  between  1.1  jum 
and  5.7  /an,  temperatures  as  low  as  500  K  could  be 
measured.  Broadband  TFP  often  suffers  from  a 
strong  flame  emission  signal.  Such  a  signal  can  be 
properly  subtracted  for  laminar  flames  [21,22], 
However,  for  turbulent  flames,  accurate  correction 
is  impossible. 

Motivated  by  this,  Bedat  et  al.  [24]  developed  a 
narrowband  TFP  technique,  which  utilizes  a  spectral 
window  centered  at  3.9  /an  to  avoid  flame  emission 


from  gaseous  molecules.  With  liquid  nitrogen  cooled 
InSb,  they  could  also  measure  temperatures  as  low 
as  500  K.  A  heated  tungsten  filament  combined  with 
a  blackbody  was  used  for  calibration.  This  approach 
is  adequate  for  non-sooty  flames.  Therefore,  the 
present  measurements  utilized  a  similar  approach 
with  the  narrowband  filter  of  Ref.  [24]  replaced  with 
a  monochromator.  In  addition,  calibration  of  the  sys¬ 
tem  was  performed  using  a  flat  laminar  premixed 
flame  and  emission  spectroscopy  as  the  primary  stan¬ 
dard. 

Experimental  Methods 

Laminar  Flat  Flames 

The  laminar,  premixed  flat  flames  were  established 
on  a  McKenna  burner.  The  premixed  methane/air 
gases  flowed  out  from  the  central  port,  which  had  a 
diameter  of  60  mm.  An  annular  coflow  of  air  (0.9  g/ 
s)  with  an  outer  diameter  of  75  mm  was  included.  A 
constant  airflow  rate  of  0.36  g/s  was  used  for  the 
central  port.  Flames  with  three  equivalence  ratios 
(< p  =  0.6,  0.54,  and  0.46)  were  established  by  vary¬ 
ing  the  fuel  mass  flow  rate.  Choked  orifice  meters 
controlled  all  flow  rates.  The  equivalence  ratio  set¬ 
tings  also  were  checked  with  sampling  and  gas  chro¬ 
matography.  The  operating  conditions  were  selected 
to  avoid  cellular  instabilities  while  providing  near- 
adiabatic  conditions.  Thermocouple  (Pt,  Pt-Rhl0% 
with  0.3  mm  bead  diameter)  traverses  were  utilized 
to  ensure  flat  temperature  profiles  and  to  measure 
the  length  of  the  homogeneous  radiating  path. 

Turbulent  Premixed  Flames 

The  turbulent  flame  burner  and  operating  condi¬ 
tions  were  identical  to  those  used  in  Refs.  [10]  and 
[25,26] .  Premixed  methane/air  flowed  through  a  15 
mm  diameter  inner  cylindrical  tube  surrounded  by 
a  102  mm  diameter  outer  airflow  casing.  The  length 
of  the  central  tube  ensured  a  fully  developed  tur¬ 
bulent  pipe  flow.  A  pilot  flame  burner  consisting  of 
90  fine  ports  with  a  diameter  of  0.2  mm  uniformly 
distributed  along  a  mean  diameter  of  20  mm  was 
mounted  at  the  exit.  A  hydrogen  (2  mg/s)  diffusion 
pilot  flame  was  used.  The  operating  conditions  in¬ 
volved  equivalence  ratios  of  0.8  and  1  for  flames 
burning  with  a  heat  release  rate  of  4.2  kW.  The  fuel 
flow  rate  was  84  mg/s  and  the  airflow  rates  were 
2000  and  1340  mg/s.  The  coflow  airflow  rate  was 
maintained  at  540  mg/s.  The  effects  of  the  hydrogen 
pilot  were  obseived  only  within  the  first  5  mm.  The 
4>  =  1,  0.8  flames  had  burner  exit  Reynolds  numbers 
of  7000  and  8700,  respectively. 

Species  Concentration  Measurements 

Measurements  of  species  concentrations  were  ob¬ 
tained  using  sampling  and  gas  chromatography  (Shi- 
madzu  GCI4A).  The  samples  were  obtained  using  a 
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quartz  microprobe  with  an  outer  diameter  of  0.5  mm 
and  an  inner  diameter  of  80  /im  using  an  evacuated 
bulb  method  [27],  The  samples  were  dried,  sepa¬ 
rated  using  Chromosorb  102  and  molecular  sieve  5A 
columns,  and  analyzed  using  a  calibrated  detector. 
The  dry-based  measurements  were  used  to  calculate 
the  equivalence  ratio  using  a  constant  C/H  ratio  and 
a  correction  for  argon  to  the  apparent  02  measure¬ 
ment.  The  uncertainties  in  the  measured  equiva¬ 
lence  ratios  were  less  than  10%. 

Measuremen  ts  of  Spectral  Radiation  Intensities 

The  spectral  radiation  leaving  diametric  paths  at  a 
fixed  height  above  the  laminar  and  the  turbulent 
flame  burners  was  collimated  by  a  308  mm  long,  4 
mm  diameter  tube  coated  with  black  paint.  A  1/8  m 
monochromator  (Oriel  model  77250)  with  adjusta¬ 
ble  slits  and  a  liquid  N2  cooled  InSb  detector  (Grase- 
by  Infrared  model  U-1933-S)  were  used.  The  output 
of  the  detector  was  amplified  and  conditioned  by  a 
preamplifier  (Graseby  Infrared  model  U-3237).  An 
optical  chopper  (Model  SR540  Stanford  Research 
Systems,  Inc.)  was  set  at  3000  Hz  for  the  turbulent 
flame  measurements.  The  output  of  the  preamplifier 
was  fed  to  a  Lock-In  amplifier  (Stanford  Research 
Inc.  model  SR510).  The  Lock-In  amplifier  was  mod¬ 
ified  to  allow  frequency  response  up  to  5000  Hz.  The 
output  signal  was  sent  to  a  laboratory  computer  via 
an  A/D  board,  CIO-DAS  16/Ml  (Computer  Board, 
Inc).  A  high  pass  optical  filter  set  at  2.5  jum  was  used. 
The  slits  were  set  to  obtain  a  bandwidth  of  AX  —  12 
nm.  The  frequency  settings  for  the  laminar  flame 
measurements  were  a  factor  of  100  lower  than  those 
summarized  above. 

A  grating  (Oriel  Model  77301)  operational  over 
the  wavelength  range  of  2.5-8  pm  was  used  to  cover 
the  spectral  range  of  interest  for  the  C02  and  HaO 
band  emission  in  the  mid-infrared  for  the  lean  pre¬ 
mixed  flames.  In  the  current  arrangement,  the  ac¬ 
curacy  and  the  precision  of  the  wavelength  adjust¬ 
ment  were  8  nm  and  1.6  nm,  respectively. 

A  blackbody  (Infrared  Industries,  Model  IR-563) 
operated  at  three  different  temperatures  (1100, 
1200,  and  1300  K)  was  used  to  calibrate  the  system, 
and  the  linearity  of  the  calibration  was  within  1%. 
Experimental  uncertainties  (95%  confidence  level) 
for  the  laminar  and  turbulent  flame  measurements 
were  less  than  2%  and  5%,  respectively.  The  mea¬ 
surements  were  repeatable  within  these  limits. 

Narrowband  Thin  Filament  Pyrometry 

A  single  SiC  fiber  with  a  diameter  between  10  and 
20  pm  was  stretched  across  a  holder  with  100  mm 
spacing  between  two  lubricated  pulleys.  Small 
weights  attached  at  each  end  of  the  fiber  kept  it  taut 
in  the  flame.  The  present  narrowband  TFP  utilized 
the  entire  spectral  radiation  intensity  measurement 


apparatus  discussed  above,  except  for  the  collimat¬ 
ing  tube.  Instead,  two  100  mm  diameter  CaF2  lenses 
were  used  to  focus  the  emission  signal  from  approx¬ 
imately  250  pm  length  of  the  SiC  fiber  onto  the  InSb 
detector.  The  monochromator  was  set  at  a  wave¬ 
length  of  X  =  3.82  pm  to  avoid  flame  emission,  with 
a  bandwidth  AX  =  0.15  pm  to  allow  adequate  signal. 

The  narrowband  TFP  has  excellent  precision  (5 
K).  Proper  calibration  of  the  fiber  and  the  associated 
optics  and  electronics  is  needed  for  obtaining  good 
accuracy.  Unfortunately,  different  researchers  have 
followed  different  calibration  techniques.  In  the 
present  work,  emission  spectroscopy  for  one  laminar 
flame  was  used  as  the  primary  calibration  method  in 
conjunction  with  a  radiation  model.  The  TFP  mea¬ 
surements  were  then  checked  against  the  adiabatic 
flame  temperatures,  which  were  calculated  using 
measured  species  concentrations.  The  net  result  was 
that  an  estimated  accuracy  of  30  K  was  achieved  with 
the  present  procedure  at  a  measurement  tempera¬ 
ture  of  1500  K  (1.5%). 

The  radiation  corrections  to  the  TFP  data  in  the 
laminar  flames  were  calibrated  out  because  of  the 
constant  convection  velocity  over  the  fiber  for  the 
different  operating  conditions.  For  the  turbulent 
flames,  an  energy  balance  over  the  fiber  involving 
heat  transfer  from  the  hot  gas  and  radiation  to  the 
black  cold  surroundings  was  performed.  An  empir¬ 
ical  correlation  between  the  Nusselt  and  the  Rey¬ 
nolds  number  was  used  [28].  The  Nusselt  number 
was  evaluated  using  the  measured  mean  velocity  and 
the  fiber  diameter  as  the  quantities  needed  for  de¬ 
fining  the  Reynolds  number.  Based  on  the  magni¬ 
tude  of  the  radiation  correction,  the  uncertainties  in 
the  correction  introduced  approximately  30  K  ad¬ 
ditional  inaccuracies  in  the  estimated  gas  tempera¬ 
tures. 

Theoretical  Methods 

We  avoided  the  uncertainties  of  the  flame  struc¬ 
ture  models  in  the  radiation  predictions  by  utilizing 
the  measured  mean  and  mis  temperatures  to  gen¬ 
erate  time  and  space  series  in  the  present  calcula¬ 
tions.  The  combined  reaction  progress  variable  and 
mixture  fraction  approach  of  Ref.  [10]  was  utilized. 

The  instantaneous  temperature  is  a  function  of  the 
reaction  progress  variable  and  the  local  mixture  frac¬ 
tion.  The  mixture  fraction  is  equal  to  unity  for  the 
material  (fuel  and  air)  flowing  from  the  burner  tube 
and  zero  for  the  coflow  air.  The  reaction  proceeds  to 
completion  prior  to  mixing  of  the  coflow  air  [10], 
This  was  applied  on  an  instantaneous  basis  here. 
Thus,  if  the  equivalence  ratio  was  equal  to  that  of 
the  incoming  mixture,  the  temperature  and  species 
mass  fractions  were  uniquely  related  to  the  respec¬ 
tive  reaction  progress  variables.  The  progress  vari¬ 
ables  based  on  temperature  and  product  species 
concentrations  were  assumed  to  be  identical. 
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=  T p  ~  T  =  YCq2  =  YhoO  (1) 

Tp  -  Tr  YC02p  yH20p 

where  c  is  the  reaction  progress  variable  calculated 
using  the  temperature  T ,  and  Yco2  and  YHo0  are  cal¬ 
culated  using  the  c  and  the  values  of  these  quantities 
in  the  product  stream  (denoted  by  subscript  P).  Tp 
is  the  peak-measured  temperature,  and  Tr  is  the  in¬ 
coming  fuel/air  mixture  temperature.  This  assump¬ 
tion  was  not  valid  in  the  flame  zone  and  for  non¬ 
unity  Lewis  numbers  because  of  the  multiple  steps 
involved  in  the  detailed  chemistry  process.  The  facts 
that  the  reaction  zones  were  relatively  thin  and  that 
the  Lewis  numbers  were  close  to  unity  for  the 
mixtures  under  consideration  may  control  the  errors 
caused  by  the  simplification. 

II  the  measured  mean  equivalence  ratio  was  less 
than  that  of  the  incoming  mixture,  then  the  species 
concentrations  were  calculated  as 


Z  = 


yH2Q 

Yh2op 


(2) 


where  Tc  is  the  temperature  of  the  coflow  air  and  Z 
is  the  mixture  fraction  defined  as  the  mass  of  the 
incoming  mixture  divided  by  the  total  mass  at  a  lo¬ 
cation  beyond  the  flame  sheet.  Use  of  detailed 
chemistry  calculations  had  little  effect  on  the  present 
conclusions. 

Once  the  species  concentrations  time  series  are 
determined  from  the  temperature  time  series,  the 
radiation  calculations  can  be  performed  if  these 
quantities  are  known  at  all  points  along  the  radiation 
path.  The  time  and  space  series  analysis  described 
the  instantaneous  temperature  fluctuation  at  a  point 
s  +  As  along  the  radiation  path  at  time  t  +  At  as 


T'(s  +  Asp  +  At)  T'(s,t  +  At) 

(T’(s  +  Asf)m  ~  <”_(T7(s)2)1/2 

+  </>2  — -  +  a(s  +  Asp  +  At)  (3) 

(T'(s  +  As)2)1'2 

The  first  two  terms  on  the  right-hand  side  represent 
the  correlated  fluctuation  in  temperature,  where  <t>i 
and  <p2  are  related  to  the  temporal  and  spatial  cor¬ 
relation  coefficients  as 


p{As)  —  p(At)p(As ,  At) 
1  -  p2(As,  At) 

p(At)  —  p(As)p(As,  At) 

1  —  pz{As,  At) 


(4) 


The  third  term  a(s  +  As,  t  +  At)  is  the  random  part 
of  the  temperature  fluctuation  whose  mean  value  is 
zero  and  the  mean  square  value  is  given  as 


az(s  +  As,  t  +  At)  =  1  —  4> i 

-  4>l  -  2(f>i<t>zp(As,  At)  (5) 

where  the  correlation  coefficients  are  given  as 


p(As) 


T'(s  +  Asp  +  At)T'{s,t  +  At) 
(T'(s  +  As)2)1'2  (T'(s)z)1/z 


p(At) 


T'(s  +  Asp  +  At)T'(s  +  Asp) 
T’(s  +  As)2 


(6) 


p(As,  At) 


T'(s  +  As,t)T'(sp  +  At) 

(■ T'(s  +  As)2)1'2  (T'(s)2)1/2 


In  the  present  simulations  the  correlation  coeffi¬ 
cients  were  modeled  as 


p(As)  =  exp(^  —  —  j,p(At)  =  exp 

p{As,At)  =  p(As)p(At)  (7) 

where  l\  is  integral  length  scale  of  turbulent  tem¬ 
perature  fluctuations  and  T\  is  the  integral  time  scale 
of  turbulent  temperature  fluctuations.  There  is  some 
support  for  the  exponential  nature  of  the  spatial  and 
temporal  correlation  coefficients  from  non-reacting 
jet  and  non-premixed  flame  literature.  However,  the 
present  extension  to  premixed  flame  radiation  is  a 
first  step,  justified  only  by  the  lack  of  relevant  infor¬ 
mation. 

Spectral  radiation  intensities  were  found  by  solv¬ 
ing  the  equation  of  radiative  transfer  covering  2.5  to 
5  p in  wavelength  range.  This  range  covers  the  two 
important  bands  of  C02  and  one  important  band  of 
H20.  Narrowband  analyses  in  RADCAL  [29]  were 
utilized.  A  more  recent  model  [30]  was  also  consid¬ 
ered  but  did  not  produce  any  improvements  for  the 
present  conditions. 

For  a  homogeneous  path  of  known  length  and 
known  components  of  the  participating  medium  in 
the  laminar  TFP  calibration  flame,  RADCAL  is  used 
in  an  inverse  iterative  procedure  to  find  the  tem¬ 
perature  and  concentrations  from  the  measurements 
of  spectral  radiation  intensities.  From  the  calculated 
C02  and  H20  concentrations,  the  equivalence  ratio 
is  also  estimated. 


Results  and  Discussion 


Laminar  Flames 

Figure  1  shows  the  measurements  and  predictions 
of  spectral  radiation  intensities  for  the  laminar  flat 
flames.  The  <fi  =  0.54  flame  was  used  for  TFP  cal- 
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Fig.  1.  Measurements  and  predictions  of  spectral  radi¬ 
ation  intensities  in  laminar  premixed  flat  flames. 


TABLE  1 

Comparison  of  measurements  and  predictions  of 
temperature  and  equivalence  ratio  for  laminar  flat  flames 


</>GC 

0.46 

0.54 

0.60 

<l>  ES 

0.48 

0.53 

0.59 

Ttit  (K) 

1384 

1498 

1591 

Tes  (K) 

1395 

1498 

1558 

Fade  (K) 

1402 

1499 

1600 

ibration.  Therefore,  the  agreement  between  mea¬ 
surements  and  predictions  for  this  flame  illustrates 
the  efficacy  of  the  inverse  procedure.  The  TFP  mea¬ 
sured  temperatures  for  the  other  two  flames  were 
used  in  conjunction  with  RADCAL  to  predict  the 
spectral  radiation  intensities.  The  agreement  (cu¬ 
mulative  error  of  less  than  1%  over  all  wavelengths) 
between  the  measurements  and  the  predictions  for 
these  two  flames  supports  the  RADCAL  model.  A 
comparison  of  the  equivalence  ratios  and  tempera¬ 
tures  is  shown  in  Table  1.  The  agreement  between 
the  emission  spectroscopy  and  TFP  data  and  the  adi¬ 
abatic  temperature  estimates  shows  excellent  inter¬ 
nal  consistency. 


Fig.  2.  Measurements  of  mean  and  rms  temperature 
distributions  at  x/D  =  2.67,  4,  and  6.67  in  a  (/)  =  0-8.  Q 
=  4.2  kW  lean  premixed  turbulent  flame. 


Turbulent  Flames 

Measurements  for  the  Q  =  4.2  kW,  (f>  =  0.8  flame 
are  considered  as  an  example.  Fig.  2  shows  the  mea¬ 
surements  of  mean  and  rms  temperatures  along  di¬ 
ametric  paths  within  the  flame  at  three  heights.  The 
peak  mean  temperatures  are  close  to  1900  K  at  all 
three  positions.  The  x/D  =  6.67  location  is  very 
close  to  the  flame  tip.  The  rms  temperatures  near 
the  peak  temperature  locations  are  veiy  low,  but  they 
increase  significantly  near  the  edge  and  near  the 
flame  tip. 

The  measurements  of  mean  equivalence  ratio  are 
shown  in  Fig.  3.  The  data  show  that  the  equivalence 
ratio  remains  close  to  the  0.8  value  for  a  wide  region 
near  the  axis  and  then  decreases  to  the  value  of  zero 
in  the  coflow.  In  the  present  model,  these  data  are 
used  to  select  either  the  progress  variable  based 
(equation  1)  or  the  mixture  fraction  based  (equation 
2)  relationship  between  temperature  and  species 
concentrations. 

Figure  4  shows  measurements  and  predictions  of 
mean  spectral  radiation  intensity  for  the  turbulent 
flame  at  the  three  locations.  The  baseline  integral 
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Fig.  3.  Measurements  of  mean  equivalence  ratio  distri¬ 
butions  at  x/D  =  2.67,  4,  and  6.67  in  a  <f>  =  0.8,  Q  =  4.2 
kW  lean  premixed  turbulent  flame. 


length  scale  for  turbulent  fluctuations  was  selected 
to  be  4  mm  based  on  past  estimates  [25,26],  The 
timescale  was  calculated  based  on  the  prescribed 
length  scale  and  mean  velocity  data  from  Ref.  [25] 
in  conjunction  with  Taylor's  hypothesis.  The  pre¬ 
dicted  mean  spectral  intensity  is  independent  of  the 
choice  of  the  length  scale  as  shown  by  the  calcula¬ 
tions  with  two  additional  length  scales  for  the  x/D 
=  2.67  and  x/D  =  6.67  locations. 

Figure  5  shows  the  rms  data  and  predictions  of 
the  spectral  radiation  intensities.  The  agreement  be¬ 
tween  the  measurements  and  predictions  with  the  4 
mm  integral  length  scale  are  excellent  for  the  x/D 
=  4  position.  However,  the  rms  is  overpredicted  at 
x/D  =  2.67  and  underpredicted  at  x/D  =  6.67.  If 
the  integral  length  scale  is  reduced  to  2  mm  at  x/D 
=  2.67  and  increased  to  8.5  mm  at  x/D  =  6.67, 
excellent  agreement  between  measured  and  pre¬ 
dicted  spectra  is  achieved.  The  monotonic  increase 
of  integral  length  scale  from  2  mm  at  x/D  =  2.67  to 


Fig.  4.  Measurements  and  predictions  of  mean  spectral 
radiation  intensity  at  x/D  =  2.67,  4,  and  6.67  in  a  <j>  =  0.8, 
Q  =  4.2  kW  lean  premixed  turbulent  flame. 


8.5  mm  at  x/D  =  6.67  is  qualitatively  consistent  with 
similarity  scaling.  However,  integral  length  scale  data 
are  necessary  to  independently  confirm  these  find¬ 
ings. 

Figure  6  shows  measurements  and  predictions  of 
PDFs  of  spectral  radiation  intensities.  Gaussian 
curves  are  included  for  reference.  The  analysis 
misses  the  highest  intensity  excursions,  suggesting 
that  the  correlated  occurrence  of  highest  tempera¬ 
tures  (flames)  is  not  captured.  However,  the  level  of 
agreement  is  quite  satisfactory  considering  that  this 
is  the  first  attempt  at  predicting  radiation  statistics 
of  turbulent  premixed  flames. 

Conclusions 

1.  Measurements  and  predictions  of  spectral  radia¬ 
tion  intensities  in  laminar  flames  show  that  the 
RADCAL  algorithm  and  constants  yield  excellent 
results  for  the  present  conditions.  Based  on  this, 
emission  spectroscopy  can  be  used  as  an  absolute 
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Fig.  5.  Measurements  and  predictions  of  rms  of  spectral 
radiation  intensity  at  x/D  =  2.67,  4,  and  6.67  in  a  r/>  =  0.8, 
Q  =  4.2  kW  lean  premixed  turbulent  flame. 


r Data  - Gaussian  —  Predictions 

1.0  . ini i| i in . . . . . . . . .  1.0 

0.8  ■  ■  0.8 


t  x/D  =  6.67  L,  =  8.5  mm  1 

0.6  I  H  0.6 


-4-3-2-10123  -3-2-10123 

g_  1  *  'mean 
‘mean 


Fig.  6.  Measurements  and  predictions  of  PDFs  of  spec¬ 
tral  radiation  intensities  at  x/D  =  2.67,  4,  and  6.67  in  a  <f> 
=  0.8,  Q  =  4.2  kW  lean  premixed  turbulent  flame. 
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COMMENTS 


Friedrich  Dinkelacker,  University  of  Erlangen,  Ger¬ 
many.  What  is  a  typical  amount  of  temperature  reduction 
due  to  radiation,  for  example,  for  a  laminar  or  slightly  tur¬ 
bulent  lean  premixed  methane/air  flame  with  <f>  =  0.60. 
Also,  what  is  the  maximum  flame  temperature  compared 
to  the  adiabatic  flame  temperature? 

Author’s  Reply.  For  a  laminar  premlxed  methane/airflat 
flame  with  (ft  =  0.6,  the  maximum  measured  temperature 
is  approximately  20  K  lower  than  the  adiabatic  flame  tem¬ 
perature.  Detailed  chemistry  calculations  considering  the 
radiation  heat  loss  show  that  radiation  from  gaseous  mol¬ 
ecules  is  the  direct  and  indirect  cause  of  this  temperature 
reduction. 


Robert  W.  Pitz,  Vanderbilt  University,  USA.  Rather  than 
modeling  the  turbulent  temperature  correlations,  you 
could  image  a  longer  silicon  carbide  fiber  and  measure  the 
instantaneous  temperature  distribution  along  the  line  of 
sight  directly.  Do  you  have  any  plans  to  measure  these  in¬ 
stantaneous  temperature  distributions  and  compare  the 
measured  correlations  to  your  modeling  results? 

Author’s  Reply.  We  are  in  the  process  of  making  such 
measurements  in  our  laboratory  now. 


